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ABSTRACT 



The acceleration cites of ultra-high-energy (UHE) protons can be traced by the footprint left by these particles propagating through 
cosmic microwave background (CMB) radiation. Secondary electrons produced in extended region of several tens of Mpc emit their 
energy via synchrotron radiation predominantly in the initial direction of the parent protons. It forms a non-variable and compact 
(almost point-like) source of high energy gamma rays. The importance of this effect is increased for cosmologically distant objects; 
because of severe energy losses, UHE protons cannot achieve us even in the case of extremely weak intergalactic magnetic fields. 
Moreover, at high redshifts the energy conversion from protons to secondary particles becomes significantly more effective due to the 
denser and more energetic CMB in the past. This increases the chances of UHE cosmic rays to be traced by the secondary synchrotron 
gamma radiation. We discuss the energy budget and the redshift dependence of the efficiency of energy transfer from UHE protons to 
synchrotron radiation. The angular and spectral distributions of radiation in the gamma- and X-ray energy bands are calculated and 
discussed in the context of their detectability by Fermi LAT and Chandra observatories. 
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1. Introduction 

Although there is a little doubt that UHE cosmic rays, E > 
lO'^eV, are produced in extragalactic sources, the nature and ori- 
gin of these objects remain highly unknown. In the case of very 
weak intergalactic magnetic fields (IGMF), the initial directions 
of these particles is only moderately entangled, and thus the in- 
formation about their acceleration sites can be partly preserved. 
The quantitative treatment of propagation of UHE cosmic rays 
1 is limited because of large (orders of magnitude) uncertainties 
■ of relevant model parameters. Especially it concerns the IGMF. 

The cuiTent measurements give only upper limit B ~ 10"^ G 
. for the component p arallel to the line of sight (Rvu et al, 1998; 
' iBlasi & OIintdrT999h^ the lower bound is estimated as small as 
(seelll 



of B ~ 10"'^ G (see^ Tavlor et al.l (1201 Ih and references therein) 
from the non-observations of GeV gamma-ray counterparts from 
distant TeV blazars. The large scale structure of IGMF is also 
very uncertai n; the attempts to model usually lea d to quite dif- 
ferent results (iDolag et al.ll2005l: ISid et al.ll2004l) . These uncer- 
tainties significantly limit the potential of the so-called cosmic 
ray astronomy. Moreover, even in the case of very weak mag- 
netic fields, the cosmic ray "horizon" of UHE cosmic ray sources 
through direct detection of these particles is limited by interac- 
tions with CMBR. For highest energy particles, E ~ 10^°eV, it 
does not exceed 100 Mpc. At lower energies cosmic rays can 
propagate to significantly larger distances, but they are dramat- 
ically deflected by galactic magnetic fields. These two factors 
significantly limit the potential of the so-called proton astron- 
omy. 

An important information about the acceleration sites of 
UHE protons is contained in gamma rays produced during prop- 
agation of protons through the CMBR. In the case of very low 
intergalactic magnetic field, the secondary products of py inter- 



actions initiate electromagnetic cascades in CMBR. The effec- 
tive development of such cascades is determined by the domi- 
nance of Compton energy losses over the synchrotron losses of 
secondary electrons. Generally, this condition is quite relaxed; it 
requires IGMF weaker than 10 '°G. However, the second condi- 
tion of detection of cascade gamma rays in the direction of par- 
ent UHE protons requires much weaker magnetic fields, smaller 
than 10" '^G. Otherwise the relatively low energy electrons, and 
consequently the GeV and TeV gamma-ray emission, will be 
immediately isotropised. Formally, such week intergalactic mag- 
netic fields cannot be excluded. Moreover, recently it has been 
argued that very high energy tails of gamma-ray spectra of dis- 
tant TeV blazars might be contributed by cascades initiated by 
UHE protons in intergalactic med ium with magnetic field as 
small as 10 '^G dEssev et al.ll201 Ih . 

Remarkably, even in the case of strong IGMF, we might 
expect physically extended (but looking point-like) gamma-ray 
sources formed around production sites of UHE cosmi c rays 
(Gabici & Aharonian 2005; Aharonian et al. 2010; Kotera et al.l 



1201 Ih . In this case gamma rays are produced through syn- 
chrotron radiation of secondary electrons from photomeson in- 
teractions. Protons lose significant fraction of their energy on 
scales of several tens of Mpc. If the magnetic field in these re- 
gions located at relatively small redshifts, z «; 1, are signifi- 
cantly large, B > lO^^G (so the energy of secondary electrons is 
released in the gamma-ray band), but smaller than lO^^G (thus 
the initial directions of electrons are not dramatically changed), 
the resulting gamma-ray synchrotron sources can be quite com- 
pact with an angular size as small as 0.1°. 

For cosmologically distant sources embedded in denser 
(wph oc (1 H- z)'^) and more energetic (T cc I + z) CMBR, pro- 
tons lose their energy on distances significantly smaller than 
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100 Mpc. In the presence of magnetic field of strength B — 
10 ''-10 **Gthe high energy electrons produced in photomeson 
processes intensively emit synchrotron radiation in GeV range 
of energies. The electrons radiate most of their energy on al- 
most rectilinear part of the path that gives, along with the small 
deflection of protons, collimated beam of gamma rays pointing 
to the acceleration site of UHE protons. Note that the consider- 
able contribution to the high energy electrons is provided by the 
gamma rays produced in decay of mesons. The interaction of 
these gamma rays with cosmic radio background radiation oc- 
curs in the regime when the most of the energy goes to one of 
the component of the created electron-positron pair 

The probability of detection of cosmic ray sources via syn- 
chrotron radiation of secondary electrons strongly depends on 
the maximum energy of accelerated protons. In the case of 
sources with redshift z <k 1, the interaction of protons with 
CMBR is eff'ective only when the proton spectrum extends to 
lO^'^eV and beyond. Since acceleration of protons to such high 
energies can be realized only at a unique combination of a few 
principal parameters, in particular the linear size of the source, 
the strength of t he magnetic field, and t he Lorentz factor of the 
the bulk motion dAharonian et al.l ( |2002|) ). the number of sources 
of 10^*^ eV cosmic rays can be quite limited compared to sources 
accelerating particles to fimax ~ lO'^eV. Since protons at high 
redshifts interact with denser and more energetic photons of 
CMBR, the requirement to /Smax is significantly relaxed, thus 
one should expect dramatic increase of the number of UHE cos- 
mic ray sources surrounded by gamma-ray halos. Furthermore, 
at high redshifts the interaction of protons with CMBR via pair 
production (Berhe-Heitler) process begins to play an important 
role and a considerable part of proton energy is converted to 
less energetic electrons compared to electrons generated in pho- 
tomesons reactions. Appearance of dense halos of Bethe-Heitler 
electrons around the source at presence of magnetic field of the 
order of B ~ 10"^ G (comparable to the field of clusters of galax- 
ies) results in radiation dominated by synchrotron X-rays. 

In this paper we study the energy and angular distribution of 
the synchrotron gamma-ray emission from cosmologically dis- 
tant sources. The calculations are based on the approach devel- 
oped in our previous work (Aharonian et al. (2010)) combined 
with cosmological eff'ects on propagation of gamma rays and 
protons. The space in the vicinity of the source, where all rel- 
evant processes occur, can be considered as conventional one. 
Therefore, taking into account denser ambient radiation at high 
redshifts, we can apply the developed formalism to calculation 
of the distribution function of gamma rays in the vicinity of the 
source. Using this function, we can easily obtain the distribution 
function of the observed radiation which has propagated cosmo- 
logical distance (see AppendixlAli. 



2. Energy budget 

An accelerator of UHE protons located at high redshift has ap- 
preciably different conditions of ambient media as compared to 
the ones in the nearby Universe. The photon field of CMB is 
denser and more energetic due to cosmological expansion. The 
increase of density by factor of ( 1 -H z)^ leads to more frequent 
interactions of UHE protons with CMB that intensifies the en- 
ergy loss rate. The average energy of photons is also increased 
by factor of (1 + z) that decreases the threshold energy of the 
interactions for protons. The energy loss rate of protons due to 
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Fig. 1. The cooling length of protons in the intergalactic medium 
due to interactions with photons of CMBR at different redshift. 
The ratio of cooling lengths at z = and z = 3 are shown in the 
inset. 



interaction with CMB at the epoch of redshift z is expressed in 
terms of the present loss rate piE) - as 

PiE,z)^(l+zfmi+z)E). (1) 

It is convenient to present this relation in terms of the cooling 
length which is shown in Fig. [T] The cubic dependence on red- 
shift leads to considerable decrease of the cooling length. In par- 
ticular, for the constant losses at highest energies, E > 10^'eV, 
it is reduced from 15 Mpc at the present epoch to ^ 0.2 Mpc 
at z = 3. The plateau of constant losses itself extends to lower 
energies due to the energy shift. At lower energies, the combina- 
tion of the eff'ects related to the energy shift and the increase of 
density jointly results in reduction of cooling length by a factor 
larger than (1 -i- z)^. Indeed, as it is seen from the inset in Fig.[T] 
the reduction of the proton cooling length can be an order of 
magnitude larger The peaks show domains where energy loss at 
z = 3 is the most intensive relative to the case of present epoch. 
If the energy cutoff in the initial proton spectrum falls into this 
domain, the advantage of energy extraction at cosmological dis- 
tances becomes more significant. Moreover, at £ = lO'^eV, the 
cooling length is reduced from > 10'* Mpc at z = to tens of 
Mpc at z = 3. It is interesting to note that the energy loss rate of 
protons of energy E < 10'^ eV at z = 3 is comparable to the loss 
rate of £ > 10^°eV protons at the present epoch. 

Fig. |2] describes the evolution of proton energy losses and 
the efficiency of their conversion to the electron component of 
secondaries with the distance to the source at different redshifts. 
The electron component includes the electrons produced through 
all channels under consideration: pair production by protons, de- 
cay of charged mesons produced in photomeson processes and 
pair production by gamma rays produced from decay of neu- 
tral mesons. As mentioned above, in the latter process almost 
all energy of the photon goes to the energy of one of the elec- 
trons. Therefore the gamma rays can be treated as electrons. 
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r, Mpc r, Mpc 

Fig. 2. The fraction of the initial energy of protons with E > lO'^ eV lost (thick lines) and converted to the energy of electrons 
(thin lines) at the distance r from the source. The injected proton spectrum is assumed power-law with an exponential cutoff 
Jp{E) = JoE-^exp{-E/Eo), with £o = 3 ■ 10^° eV. 



Then the energy taken away by neutrinos is the difference be- 
tween the energy lost by protons and the energy converted to 
electrons. As can be seen from Fig.|2] the protons with initial en- 
ergies E > 10'^ lose only 18% of their total energy after passing 
r = 100 Mpc at the redshift z = 0.1, whereas at redshift z = 3 the 
same protons lose 93% of the available energy already at r = 70 
Mpc, where the saturation begins. It can be explained by the fact 
that the cooling length of protons in the range E > 10^^ eV do 
not exceed tens of Mpc at z = 3, while for redshifts z <s 1, the 
protons have the cooling length of cosmological scales relative 
to pair production process. As the contribution of pair produc- 
tion process increases, the share of energy lost by protons that 
goes to electrons increases from 74% at z = 0.1 to 86% at z = 3 

The fraction of the initial energy of the protons with E > 
lO'^ eV that can be converted to the electrons generated at pho- 
tomeson processes depends strongly on the position of the cut- 
off energy. Left panel of Fig. |3] shows that the efficiency of ex- 
traction of proton energy and its transfer to this component of 
electrons grows with the redshift, and at z = 3 all available for 
conversion energy is transferred at first 5 Mpc. The fraction of 
proton energy taken away by neutrinos in photomeson processes 
is 42% independently of redshift. Due to more energetic photons 
of CMBR, the threshold of photomeson interactions is shifted to 
lower energies that leads to the increase of extracted energy from 
6% to 17% (see Fig. |3]l. However, as pair production begins to 
play a significant role at high redshifts the share of the electrons 
generated in photomeson processes is decreased, from 49 % to 
22 % at the distance r = 100 Mpc. 

The importance of the shift of the reaction threshold on 
production of electrons becomes more obvious if the protons 
with energies close to threshold and higher are taken into ac- 
count. The percentage of initial energy of protons with energies 
E > lO'** eV which is converted to the energy of electrons is 
presented in the right panel of Fig. |3] While at z = 0.1 the pro- 
duction of electrons due to photomeson processes dominate at 
all distances smaller than 100 Mpc, at high redshifts the pair 
production prevails. 



3. Gamma-ray source 

Having the highest energy among secondaries, the electrons gen- 
erated in photomeson processes emit almost all their energy 
via synchrotron radiation. Therefore, because of the shift of 
the threshold of photomeson interactions, the intensity of syn- 
chrotron gamma rays is increased with the redshift of the source. 
On the other hand, the reduction of cooling length of protons re- 
sults in reduction of the apparent angular size of the region emit- 
ting synchrotron radiation, in addition to the diminution because 
of geometrical factor. 

Fig. |4] shows spectral energy distribution of synchrotron 
gamma rays observed within different polar angle in the di- 
rection to the source of protons. The geometry of expanding 
Universe leads to more sophisticated dependence of apparent an- 
gular size on distance compared to the ~ l/r dependence for the 
local Universe. It is useful to define a reference angle 61 which 
is equal to the angular size of the region with the radius r = 0.5 
Mpc located at the redshift z. Expressing the angles in units of 
corresponding reference angles, we can compare the regions lo- 
cated at different redshifts eliminating the geometrical factor. As 
expected, at large redshifts the reference angle increases with z 
(see AppendixlAl Fig. lA. 11 1. In Fig.|4]for each redshift the fluxes 
are presented within three polar angles which differ from each 
other by factor of 3 and are expressed in the units of correspond- 
ing reference angle. The maximal angle indicated in plots is the 
polar angle within which the total flux is observed. Comparison 
of the maximal angles in units of reference angles shows a ten- 
dency of decreasing the angular size of the region of secondary 
synchrotron radiation with redshift, from 90* j to 6^9. 

The interaction of synchrotron gamma rays with extragalac- 
tic background light (EBL) leads to considerable attenuation of 
the flux. As it is seen from Fig.|4]the intergalactic absorption be- 
comes substantial, depending on the distance to the source, from 
TeV energies down to tens of GeV energy range. For calculation 
of the absorbtion of gamma rays, the model of EBL developed 
in lFranceschini et al.l (l2008h have been applied. 

At high redshifts, the electrons from photomeson reactions 
are produced close to the acceleration cites of protons. In this 
case a considerable part of emitting electrons might be found 
in much stronger magnetic field compared to the average in- 
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Fig. 3. Left panel: the fraction of the initial energy of protons with energy E > 10'^ eV converted to the energy of electrons from 
photomeson production at distance r from the source. Right panel: the percentage of the initial energy of protons with energy 
E > lO'^ eV converted to the energy of electrons from photomeson production (thick lines) and electrons from pair production (thin 
lines) at the distance r from the source. 
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Fig. 4. Spectral energy distribution of gamma rays observed within different angles in the direction to the source with EBL ab- 
sorbtion (thick lines) and without one (dashed lines) indicated for maximum angles. 0! is apparent angular size of the region with 
radius 0.5 Mpc from the distance with redshift z- The angles are specified in the units of corresponding ffl and differ from each other 
by factor of 3. The maximum angles for every redshift, the total power of injection of protons L and the beaming angle Oheam are 
indicated. 



tergalactic field. The energy spectrum of synchortron radiation 
of secondary electrons is shifted linearly with change of the 
strength of magn etic field to keep the ratio EyjB constant (see 
lAharonian et al I (120 lOt) ). Therefore the increase of the strength 
of magnetic field leads to the shift of radiation spectrum towards 
higher energies; if the radiation extends to TeV energies, the in- 
tergalactic absorption becomes quite violent; this results in the 
dissipation of almost the entire radiation. The absorbed gamma 
rays initiate cascades which contribute to the diffuse extragalac- 
tic background radiation. 

Keeping in mind the sensitivity of instruments such as 
Fermi, the detection of the collimated synchrotron radiation 
from sources located at high redshifts is possible only in case of 
very powerful AGN. The anisotropic injection of UHE cosmic 



rays allows to reduce the required power of source. Until deflec- 
tion of protons is smaller than the angle of a jet, there is no differ- 
ence between the spherically symmetric and anisotropic case. In 
both cases the observer will see identical pictures. For the given 
power of injection the existence of anisotropy with the opening 
angle 9beam of jet means the increase of the flux of gamma rays by 
the factor 4/0^^^^^^ compared to the spherically symmetric case. 
The fluxes of synchrotron gamma rays in Fig. |4]are calculated 
for the case of 9beam - 10°. The deflection of protons is smaller 
than this value of the opening angle while protons produce sec- 
ondary particles. It is seen from Fig.|4]that at the power of proton 
injection L - 10'^^ - IC'^erg/s, the fluxes of gamma rays can ex- 
ceed the sensitivity (minimum detectable fluxes) of Fermi LAT 
(lAtwood et al.ll20()9l) The resolution of the Fermi LAT varies in 
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Fig. 5. Energy flux distribution of gamma- and X-rays observed 
within different angles in the direction of the source for the case 
of homogeneous magnetic field B = 1 nG. Other parameters are 
the same as in previous figures. 



dependence of photon energy from 42' at £ = 10'' eV to 4.2' 
at £ = 10" eV dAtwood et al.ll2009h . Already for the source lo- 
cated at z = 0.2 the angular size is 46', thus the source of gamma 
rays would be seen as extended one only if z < 0.2. 



4. X-ray emission 

Protons lose part of their energy via pair production. The pair- 
produced electrons have lower energies compared to the elec- 
trons produced in photomeson processes and emit synchrotron 
radiation in the lower band of spectrum. At high redshifts the 
mean free path of protons relative to pair production process is 
decreased and generated electrons are located in more compact 
region. Moreover, the fraction of energy of protons converted to 
the secondary electron component is increased and reaches to 
68% atz = 3 (see Fig. 121 IS. 

Synchrotron radiation for the strength of IGMF of 
B = 10 G and inverse Compton (IC) scattering give equal con- 
tribution to electron losses at energy E ^ 3 ■ 10'^ eV. At lower 
energies electrons lose their energy predominantly through IC 
scattering. Scattering is carried out at Klein-Nishina regime and 
almost all energy is transferred to photon. In turn the high energy 
photon produce electron-positron pair due to interaction with 
CMBR and most of the energy goes to one component of the 
pair which again suffer IC scattering. This process can be con- 
sidered as alternation of the particle state with gradual reduction 
of energy. The electromagnetic cascade produces a huge halo of 
gamma rays with energy in GeV region. Region of synchrotron 
radiation is more compact, and electrons radiate it at the place of 
their generation. In spite of this the region is still quite extended 
as can be seen from angular distribution on Fig. |5] The flux of 
synchrotron radiation drops at low energies and becomes very 
small in X-ray region. The calculation of fluxes presented on 
Fig. Intakes into account only homogeneous IGMF with strength 
B - lO^^G. However, the strength of magnetic fields close to 
the accelerator can be much higher To take this into account, we 
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consider the following spatial distribution for IGMF magnetic 
field: 

B = +Bo, (2) 

where B^i = lO^^G is magnetic field in the cluster of galaxies 
with characteristic size of ro = 1 Mpc, Bq = 10"^ G is IGMF. 
Such a dependence of the magnetic field has been chosen to have 
a dipole behavior at large distances. For this case the energy and 
angular distributions of gamma- and X-rays are presented on left 
panels of Fig. |6] and Q This magnetic field gives considerable 
increase of the fluxes of synchrotron radiation at low energies 
which is generated at small region close to the source and, there- 
fore, has narrow angular distribution. Right panels shows that 
the contribution of electrons produced in pair production pro- 
cess (lower lines) is dominant in X-ray region, whereas gamma 
rays are generated predominantly by electrons produced in pho- 
tomeson processes. Horizonta l segments indicate se nsitivity of 
Chandra dLehmer et al.l 120051: iRomano et al.l l2008l) for corre- 
sponding maximum angles of observation. As can be seen the 
low energy part of X-rays (0.5-2 keV) is detectable for specified 
power of injection of protons, whereas the radiation at higher en- 
ergies can be missing for sources located at high redshifts since 
Chandra is less sensitive in the range of energies 2-8 keV. 

As follows from Fig.|6]and|2l the electrons produced in pho- 
tomeson processes provide the main contribution to synchrotron 
gamma rays, although the most of energy lost by protons is con- 
tained in the low energy component of secondary electrons from 
the pair production process. The radiation of this component can 
be detected only in the case of a quite large, Bd - 10"^ G mag- 
netic field around the source. This can happen if the source is 
located inside a cluster of galaxies with a typical size of 1 Mpc. 
Nevertheless, even in the case of intensive pair production that 
takes place in cosmologically distant objects, only a small frac- 
tion of pairs is produced in the relatively compact region with a 
linear size of ~ 1 Mpc (see Fig.|2|; the most of energy of protons 
is converted to the energy of extended gamma ray halo. 

5. Summary 

High energy gamma rays are unavoidably formed around the 
sources of UHE cosmic rays because of synchrotron radiation 
of secondary electrons produced at interactions of highest en- 
ergy protons with the cosmic microwave back ground radiation. 
In our previous paper (lAharonian et al.l (1201 Ob ) we have shown 
that even for relatively large intergalactic magnetic fields in the 
neighborhood of UHE cosmic ray accelerators, B ~ 10"^ - 
10"^ G, these process lead to formation of high-energy point- 
like gamma-ray sources. Because of relatively weak gamma ray 
signals, the chances of detection of such sources are higher for 
objects located in the nearby Universe, namely at distances less 
than 100 Mpc. Since the efficiency of conversion of the pro- 
ton energy to secondary gamma rays is dramatically reduced at 
protons energies E < 10^" eV, one may hope to detect gamma 
rays only from extreme objects accelerating protons to energies 
10^" eV and beyond. Given the fact that this requirement can be 
satisfied only in the case of unique combination of parameters, 
as well as the limited volume of the < 100 Mpc region, realis- 
tically one can expect not a very large number of such sources. 
One can gain a lot if extends the study to cosmological distances. 
First of all, this allows to probe the most powerful objects in 
the Universe (e.g. quasars and AGNs) in which more favorable 
conditions can be formed for acceleration of protons to energies 



6 



A.Yu. Prosekin et al.: Non- variable cosmologically distant gamma-ray emitters 



10 



12 



\ 10-^3 
o 

LU 



10 



16 



. z=0.1 


1=10^*^ erg/s 




^ z=0.5 




: 










^/^^ - 


: ^ 























9 /_9-|4' : 

"0.1 vj. 1 1 






e„5*=2.74' : 


3^ 10^ 10"^ 


10^ 10^ 


10^ 10^ 




E, eV 






Fig. 6. Energy flux distribution of gamma- and X-rays observed within different angles in the direction of the source (left panel) and 
within maximum angle (right panel) for the case of two-band magnetic field. For the right panel, the total radiation from photomeson 
and pair production electrons (upper lines) and the radiation from pair production electrons (lower lines) are indicated. Horizontal 
segments present the Chandra sensitivity for the corresponding maximum angles of observation. The other parameters are the same 
as in previous figures. 




10 eV. Secondly, because of higher temperature of the CMBR 
at cosmological epochs, i.e. because of denser and hotter relic 
photons, less energetic protons (with energy down to 10^' eV) 
can lead to effective production of gamma rays. Given that the 
conditions of acceleration of protons to 10'^ eV in suspected 
cosmic accelerators are muc h relaxed as compared t o the 10^*^ eV 
extreme accelerators (see Aharonian et al.l (120021) '). we should 
expect dramatic increase of such gamma-ray sources. Another 
factor of enhancement of number of sources comes from the in- 
crease of the volume of the explored region to redshifts z > 3. 
An obvious caveat in this case is the decrease of gamma ray flux. 
However, this factor can be at least partly compensated by the 
huge power of objects in the remote Universe. Moreover, due to 
dramatic reduction of mean free paths of UHE protons at cosmo- 
logical epochs z > 1, the conversion efficiency of proton energy 
to gamma rays is increased almost an order of magnitude, which 
makes these objects an extremely effective gamma-ray emitters. 



Finally, since the secondary gamma-ray emission generally fol- 
lows the direction of parent protons, the beamed cosmic acceler- 
ators like GRBs and blazars seem to be quite attractive targets for 
search of point-like but steady GeV gamma-ray emission from 
cosmologically distant objects. 

Appendix A: Angular size of sources at large 
redshifts and distribution function in the 
expanding Universe 

If the radiation is generated on scales smaller than cosmological 
ones, the relevant processes can be considered as they take place 
in the conventional stationary Euclidean space. However, when 
radiation propagates over cosmological distances, the expansion 
of space should be taken into account. Based on the law of free 
motion of massless particles in expanding space, one can find the 
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distribution function of photons at the observation point. Then 
we need to convert the distribution function calculated in the 
reference frame with origin at the so urce to the reference frame 
of the observer (iKelner et alJ (l201 lb ). In the present paper the 
free propagation of photons is considered in the flat expanding 
Universe with parameters Qa = 0.73, Q,,, - 0.27 and Ho = 71 
km/s Mpc. 
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Fig. A.l. Angular size of the source with diameter Z) = 1 Mpc located 
at different redshifts z. 



We consider an isotropic gamma-ray source of radius lo- 
cated at redshift z. Let us assume that photons escape this region 
with spherically symmetric distribution fziE, 0), where E is the 
energy of photon and 6 is the angle between photon momen- 
tum and radial direction at the escape point. Finally we assume 
that after the escape gamma rays propagate freely. In the case 
of small angles the distribution function /o(z, E, Q.(6)) of gamma 
rays at the observation point integrated over the solid angle Q 
(with polar angle ff) can be expressed in the following form: 

/o(z,£,Q(0)) = 27rJ^j j ME(l+z),x)xdx, (A.l) 



where 



0, = 



dz' 



(A.2) 



The final result (lA.lb differs from analogous one corresponding 
to the stationary space by the dependence on the redshift as well 
as by the nonlinear dependence of angular size on distance. The 
latter is determined by the parameter 0, given by Eq. ( IA.2I ). The 
parameter AO - 20, has the meaning the angular size of the 
isotropically emitting source. In the case of anisotropic source, 
the angular size of radiation cannot be arbitrary large and is lim- 
ited by AO. In Fig. lA.ll we show the angular size of the emission 
region with diameter D - IR, = 1 Mpc. The parameter AO has 
a minimum 1 .95 arcmin at redshift z = 1 .64 when reces sion ve- 
locity equals to the speed of light c (iKelner et al.l (1201 ih ). 
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